Suspicion on Engrafting HBT From Astronomy to Heavy Ion Collision 
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HBT method in astronomy and heavy ion collision is contrasted in present article. Some differences 
are found and validity of using HBT in heavy ion collision is suspected. 
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I. INTRODUCTION 



The method of two-particle intensity interferometry was discovered in the early 1950's by Hanbury Brown and 
Twiss (HBT) |l| who applied it to the measurement of the angular diameter of stars and other astronomical objects. 

Then several authors [S IE U E3 have proposed HBT studies to probe source structure in heavy ion collision. 

The spirit of this method is reviewed and explained by Ulrich Heinz |2|: "two random point sources a and b on a 
distant emitter, separated by the distance R, emit identical particles with identical energies E p — (m 2 +p 2 ) 1 ^ 2 which, 
after travelling a distance L, are measured by two detectors 1 and 2, separated by the distance d (see Figure^). L 



(a) (b) (c) 




FIG. 1: Measurement of the separation R of two sources a and b by correlating the intensities in detectors 1 and 2 at varying 
distances d. a: The general scheme, b: The specific situation in astronomy, c: The specific situation in particle physics. 

should to be much larger than R or d. The total amplitude measured at detector 1 is then 

Ai = -(a e <&T+*<0 + (3 e ^ prib+< t"'^ , (1) 

where a, (3 are the amplitudes emitted by points a and b, </>a,</>6 are their random phases, ri a ,ri& are their distances 
to detector 1. For some time ,the two-particle intensity correlation function is thus given by 



C(R, d) = = 2 = 1 + ' ' cos [p r lQ - r 2a - r lb + r 2b ) . (2 

{h){h) (A 1 Al)(A 2 Al) {\oi\ 2 + Wr ' ' 



For large L ;§> R,d, the argument of the second, oscillating term becomes 

r\a - r 2a - r u + r 2 b — > (cos(d, R) - cos(d, L) cos(i?, L)^j . (3) 

Note the symmetry of this expression in d and R, the separations of the detectors and of the emitters." this symmetry 
is lost in the two practically relevant limits: 

f. In case of R ti,the cosine-term in (|2~j) reduces to cos(ci*(p a — p^)), with Pa,b — P&a,b- 

2. In case of R <C d,Then (see Figure^;) the cosine-term in © becomes cos(i?-(pi— p 2 ))- 

So we can get the intensity correlation function 

C( P i -p 2 ) = l + J d 3 Rp(R) cos(i? • (ft - p 2 )) . (4) 

for a continuum of a static sources described by a distribution p(R) of their relative distances. 
The upper is deduction in theory.Now let's see how to get the correlation function in experiment. .3]" Experimentally, 
the two-particle correlation function is obtained from the ratio C 2 (q) = A(q)/i?(q) (normalized to unity at large q), 
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where A(q) is the measured two-pion distribution of pair momentum difference q = P2 — Pi, and B(q) is the mixed 
background distribution Q, calculated in the same way using pairs of particles taken from different events." Note 
P2,Pi is the primary momentum of particle, i.e. the momentum of the particle when the particle is nearest to the 
source. It corresponds to the momentum when the particle leaves the source in theory. 

Really the variable used in correlation is a Lorentz invariable <3mi>.I3"Qmtj = \J (pi — P2) 2 — (B\ — E2) 2 " -After 
transforming and C2(q) into function of Qinvjb is believed that relationship between experiment and theory is 
established. 



II. ESSENTIALS IN HBT AND SUSPICION 



Correlation function in @ contains five variables R,d,p±,p2 and L. Correlation function in Q contains four variables 
d,pi,p2 and L. The diminishing of R is caused by the integral over i?.The process using HBT equals the process 
retrievinng these variables from experiment. 

In astronomy, d is known quantity , L can be got from other astronomical method. Although measuring pi and p 2 is 
difficult, we can cancel them by integral over them that is recording all particles entering the detector in experiment. So 
all variables are got from experiment and HBT succeeds in astronomy. 

In heavy ion collision ,in order to retrieve the correlation function from experiment, these four variables d,pi,p2 and 
L must be fixed. We can get p\ and P2 from detector easily.The following and last question is to get d and L,in other 
words, where the position the particles are detected is. The thought that we can cancel d and L after the integral 
over them in experiment is wrong because we can't place point detectors throughout the space and even if we place 
point detectors throughout the space in ideal we still can't make this integral in experiment just because that particle 
entering the second detector after entering the first detector has changed its quantum property by the first detector 
so becomes invalid. Another opinion regarding the first point left in TPC as the position the particle be detected 
is also wrong because that the reconstructed track of particle is classical and the point is on this classical track ,as 
well as Q becomes invalid. 

If d and L can't be got from experiment, HBT can't be used in heavy ion collision. Now let's see what we got from 
experiment in Q.The process in experiment is :firstly , determining the momentum and position of particle by the 
points left in detector;secondly, reversing the particle along the track, by the way the classical concept of orbit are 
needed, to the nearest point in the track to the source and calculating the momentum at the point;thirdly,obtaining 
the correlation function C2(q) = A(q)/B(q) (normalized to unity at large q), where A(q) is the measured two-pion 
distribution of pair momentum difference q = P2 — Pi, and B(q) is the mixed background distribution calculated 
in the same way using pairs of particles taken from different events. From these process we can know the momentum 
used in calculating the correlation function is the momentum when particle is in the point that is nearest to the 
source. This point corresponds to the point in which the particle emits from the source in theory. 

Now we can reply the question how to get d and L,in other words, where the position the particles are detected is 
.In Q,the particles are detected on the surface of the source and L=0 ,d ~ 0. So the deduction in f| is disabled here. 



III. FORMULA FOR EXPERIMENT 



Now let's see what we have done in data analysis in 3]. Define P(9, (f, r) is the emitting probability of the source 
where 9, ip is polar angle and azimuth angle respectively and r is the proper time of the source. For an event with 
M particles considered, n s ,i (q,q + dq) is the number of particles whose momentum difference with the i-th particle 
locates q and q + dq, where q is q = P2 — Pi or qi nv in different i (q, q + dq) is the number of particles in 

other event whose momentum difference with the i-th particle locates q and q + dq, it is like n s ,i (q, q + dq) expect 
that n m ,i (q, q + dq) is for different events whereas n 8 ,i (q, q + dq) for same event. 

For a completely random source , such as an artifical source produced by computer obeying the distribution 
P(9, <p, t) of real data, it is easy to know 



n s ,i(q,q + dq) =n mn (q,q + dq) (5) 

C 2 (q) = gg°^4li±^ = 1 (6) 
Ei=o n m , t (q,q + dq) 

in which average is over all events. By now, it is clear that the correlation function has no direct relation with the 
emitting probability P(9, <p, t). 
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In real data, if there are some mechanisms can change n s ,i {q,q + dq) to n s ,i (q,q + dq) + n c ^ (q) , whereas the sum 
number of particles and the distribution have not changed so n m;i (q, q + dq) have not changed, we can get 

n i„\ J2tLa[ n s n (q,q + dq)+n cn (q)} 

L 4=0 n m,i (q,q + dq) 

n m ,i (g, q + dq) can be got from simulation by computer after knowing the distribution P(0, ip, r).n c ,i (q) is a variable 
connecting to physical process on the surface of the source i.e. what happens when particle is produced from partons. 
Let q = q = p 2 — Pi,0 become the result in 



IV. MECHANISMS PRODUCE n c „ (g) FOR SMALL q 

As what are interesting is the performance of correlation function at q ~ 0,i.e. the difference of momentums is 
small, only this part is discussed here. 

The meaning of n c ,j (q) is the variance of number of particles whose momentum p satisfies p—pi — q from completely 
random distribution. If the i-th particle was produced, the probability of production of particle whose momentum is 
close to pi would increase. In this case n c ,i (q) > O.Now if the i-th particle was produced, that probability would 
decrease. In this case n c ,i (<?) < 0. 

Let's see some mechanisms making n c ,i (q) positive or negative. On the surface of the source, if a fermion is 
produced, it will hold a state of specified momentum and prevent particle of this kind produced after it entering 
this state. Qualitatively,?^^ (q) will become negative. Interaction will also affect n c ,,; ((j).If there an attractive force 
between a kind of particles, n Cli (q) will become positive, vice versa, attractive force leading negative n C) , (q). 



V. CONCLUSIONS 



The most important problem in engrafting HBT from astronomy to heavy ion collision is the position where particle 
is detected can not correspond to corresponding variables in quondam theory.The momentum used in HBT is that 
when particle leave the source. At this point L—0 ,d ~ 0. So the deduction is disabled for experiment. 
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